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("■^ ' Abstract 

^N ' In [BW07| Buczyiiska and Wisniewski showed that for the Jukes Can- 

tor binary model of a 3-valent tree the Hilbert polynomial depends only on 
the number of leaves of the tree and not on its shape. We ask if this can be 
generalized to other group-based models. The Jukes Cantor binary model 
has Z2 as the underlying group. We consider the Kimura 3-parameter 
model with Z2 x Z2 as the underlying group. We show that the gener- 
alization of the statement about the Hilbert polynomials to the Kimura 

f) , 3-parameter model is not possible as the Hilbert polynomial depends on 

^^ • the shape of a 3-valent tree. 

■3 









1 Introduction 



Phylogenetic algebraic geometry studies complex algebraic varieties arising from 
evolutionary models in biology (see for example |ESSR05| . [SSOSj V The Jukes 
Cantor binary model is the simplest model. Buczynska and Wisniewski showed 
^ I in [BW07) that the Hilbert polynomial of the ideal of the Jukes Cantor binary 

model with a 3-valent tree T depends only on the number of leaves of T and 
not on the shape. We are interested if this property of the Hilbert polynomial 
^>f-\ I can be generalized to more complicated models. 

I • I The Jukes Cantor binary model is a group-based model with the underlying 

^—v I group Z2 . The most natural generalization seems to be the Kimura 3-parameter 

f^ I model, which is a group-based model with the underlying group Z2 x Z2. How- 

ever, in this paper we conclude that this generalization is not possible. For the 
Kimura 3-parameter model the Hilbert polynomial also depends on the shape of 
. a 3-valent tree. We show that the ideals of 2 different trees with 8 leaves - the 

K> I 3-caterpillar and snowflake trees (see figures on page 3) - have different Hilbert 

>^ . polynomials. The Kimura 3-parameter model being the closest model to the 

DNA binary model, it is unlikely that the property about Hilbert polynomials 
would hold for other models. 

In Section 2 we recall the construction for the Kimura 3-parameter model. In 
Section 3 we show that the Hilbert polynomials of the ideals of the 3-caterpillar 
and snowflake trees have different values when evaluated at 3 and hence the 
Hilbert polynomials are not the same. The main idea is to decompose the 
original trees to smaller trees and use toric fiber products introduced by Sullivant 
in [Sul07] . Michalek showed in [MiclOj that the Kimura 3-parameter model is 
normal, so we reduce the problem of evaluating the Hilbert polynomials of toric 
varieties to evaluating the Ehrhart polynomials of the corresponding polytopes. 



Computations are done with the help of polymake I J Gj . |JMP09j and Normaliz 

m- 

I would like to thank Christian Haase and Andreas Paffenholz for their ideas 
and technical help. 

2 Kimura 3-paranieter model 

We will not give the parametric construction of the Kimura 3-parameter model 
coming from biology here (see |SS05| ). but will directly define the toric ideal 
and the corresponding polytope of the Kimura 3-parameter model. By doing 
this we follow [SS05] and [SuIOT] . 

Let T be a 3-valent tree with n + 1 leaves labeled by 1, . . . , n + 1, let the 
root be at the leaf n + 1, and direct the edges away from the root. A leaf I is a 
descendant of an edge e if there is a directed path from e to L Denote by de(e) 
the set of all descendants of the edge e. 

For a sequence gi , . . . , 5„ in Z2 x Z2 , we define 

z(^de(e) 

where e is an edge of T and the subindices i denote simultaneously leaves and 
their labels. Let 

K[q] = K[qg,^,„^gjg^ G Z2 X Z2] and K[a] = K[a|f^le S E{T),he Z2 x Z2] 

and consider the ring homomorphism 

(j)T : K[q] -^ K[a] 

1ai.-,ar.^ n 4^ 

eeE{T) 

Definition 1. Let the ideal It — kei{(pT) be the ideal of the Kimura 3- 
parameter model with tree T. 

The ideal It is a toric ideal and we can define the corresponding polytope. 
Definition 2. Let the polytope 

Pt = conv({a e z^('^)x(^^x^^)|a" = 0Tfei,...,9J, %„...,,,„ e K[q]}) 
be the polytope of the Kimura 3-parameter model with tree T. 

Definition 3. Let the lattice Lt C z^C^^^'^^xZa) be the lattice generated by 
the vertices of Pt ■ 

Since T is an acyclic directed graph, there is an induced partial order on 
the edges of T. Namely e < e' if there is a directed path from e' to e. Let T 
be a tree that contains an interior edge e. Then e induces a decomposition of 
T as T+ * T~ where T~ is a subtree of T consisting of all edges e' £ T with 
e' < e and T^ consists of all edges e' £ T with e' ^ e. Thus r+ and T~ overlap 
in the single edge e. We root T~ by the tail of e, and keep the root of T+ 
at the original root n -I- 1. Let /^+ and /J_ denote the ideals of the Kimura 



S-parameter model with trees T+ and T~ and let K[(7]+ and K[(7]_ denote the 
ambient polynomial rings, respectively. 

Using toric fiber products of ideals (see |Sul07) ) Sullivant stated the following 
theorem. 

Theorem 4 (Sullivant |Sul07) ). Let T he a tree with an interior edge e, and 
resulting decomposition T = T^ * T~ . For each variable Qg — qgj^,...^g^ in 
IC[q],IC[q]+ and IC[q]-, let deg{qg) = Cg^. Then 

It = It+ ^ A It- 
with A = {e(o,o), 6(0,1), 6(1,0), 6(1,1)}. 

3 Counting lattice points 

1. Since the polytopes of the 3-caterpillar and snowflake trees are too large 
to compute their lattice points directly, we decompose them into smaller 
trees like shown on the figure below. 



e 
3-caterpillar 




Henceforth we use the abbreviations 3c, sn, 31, 41 for the 3-caterpillar, 
snowflake, 3-leaf and 4-leaf trees, respectively. 

In the decomposition of the 3-caterpillar tree define deg(<7g) — Cg^ for q^ 
in /C[9]3(, and A^[9]4j. Then 

^3c = ^41 ^ -4 Il\ 
with A = {6(0,0), 6(0,1), 6(1,0), 6(1,1)}. 

In the decomposition of the snowflake tree define deg(gg) = Cg^ ,g^ for 
(7g in lC[q]^^ and /C[g]4j and define deg(gg) — Cg^ with i G {1,2} for gg in 
/C[g]3j. Then 

T _ Tei,e2 y . rei re2 

-'sn — -'41 ^.A -fai ^AJ^^i 

with A = {6(0,0), 6(0,1), 6(1,0), 6(1,1)}. 



2. Denote the multigraded Hilbert function of K[(7]/J by h{K[q]/I;u). In 
[Sul07j Sullivant gives a formula for computing multigraded Hilbert func- 
tions of toric fiber products. Applying this to the decompositions of Step 
1 gives for u,w e Z^^""^^ 

h{K[q]Jh,;u) = h{K[q]JIt,;u)h{K[q]JIt,;u), 
hiK[q]JI,^;u,v)^h{K[q]JQ^'^;u,v)h{K[q]JI!^l;u)h{K[q]JIIf;v). 

3. A monomial having multidegree u E 2^2 xZa j^g^g total degree J2hez xz "'«■ 
Thus single graded Hilbert functions can be computed using multigraded 
Hilbert functions 



hiK[q]Jh,;n 
h{K[q]JI,^;n)^ 



J2 h{K[q]Jh,;u), 
J2 h{K[q]JI,n;u,v). 



4. In [MiclO] Michalek shows that the polytopes of the Kimura 3-parameter 
model polytopes are normal, hence corresponding Ehrhart and Hilbert 
polynomials are equal (see for example |Stu96j l Thus h(K[q]/lT; u) counts 
lattice points in the lattice Lt of the X!/iez xz ^h dilation of the polytope 
Pt intersected with hyperplanes {xf^ = Uh};h G Z2 x Z2. Using Step 2 
and Step 3 we get 



ehrpj^ (n) 
ehrp^^ (n) 



u.v.Yl Uh=n,J2'"h=n 



^31 



5. Using polymake and Normaliz we can compute \3PTC\{xf = u;}nLr| for 
3-leaf and 4-leaf trees. It is important to do the basis transformation be- 
fore counting lattice points, since these programs assume that the lattice 
is the standard lattice. Using formulas from Step 4 we got that in the 3rd 
dilation the 3-caterpillar polytope has 69324800 and the snowflake poly- 
tope has 69248000 lattice points. Hence their Ehrhart (and thus Hilbert) 
polynomials are different. 

Remark. Similar computations show that in the 2nd dilation the 3-caterpillar 
and snowflake polytopes have both 396928 lattice points. 
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